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Abstract. We analyze the low frequency part of charmonium spectral functions on 
Q ■ large lattices close to the continuum limit in the temperature region 1.5 < T/Tc ^ 3 as 

5 , well as for T ~ 0.75Tc. We present evidence for the existence of a transport peak above 

Tc and its absence below Tf,. The heavy quark diffusion constant is then estimated 
using the Kubo formula. As part of the calculation we also determine the temperature 
dependence of the signature for the charmonium bound state in the spectral function 
and discuss the fate of charmonium states in the hot medium. 
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P^r ! 1. Introduction 

Experimentally a substantial elliptic follow of heavy quarks has been observed at 
RHIC [1]. Various phenomenogical model studies suggest heavy quark diffusion 
^ ■ coefficient D < 1/T to accomodate data (see e.g. Ref. [2]). Theoretically the heavy 

c^ ■ quark diffusion coefficient D has been calculated by perturbative QCD in both leading 

and next-to-leading order as well as from AdS/CFT correspondence. At a^ ~ 0.2, 
leading order pQCD calculation gives 2txTD ^ 71.2 [2J while next-to- leading order 
calculation gives 2ttTD ^ 8.4 [3]. In the strong coupling limit 2itTD = 1 is obtained 
from AdS/CFT correspondence [^. Under such a circumstance a non-pertubative 
computation of heavy quark diffusion is needeqj|. 

Through Kubo formula, the heavy quark diffusion constant D relates to the vector 
spectral function as D = ^^^ lim^^^o ^^^-^ °^'^^' , where xoo is the quark number 
susceptibility and ay{u,T) is a mesonic spectral function in the vector channel. In 
the non-interacting case (Ty{uj) has a uj6{u) term [6] and consequently gives an infinity 
diffusion constant, while in the interacting case u6{uj) will be smeared into a Breit- 
Wigner form [7] and leads to a finite diffusion. The fate of quarkonia states at finite 

J Here we focus on extracting diffusion coefficient from mesonic spectral functions, other ways can be 
found in Ref. [5]. 
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temperature, which was suggested as a useful probe of QCD medium properties |S], can 
also be signaled by the deformation of spectral functions. 

The mesonic spectral function is not directly accessible through lattice QCD 
simulations and can be obtained from the inversion of the following equation 

cosh(aj(r — ^)) 



Gh{t,T) 



duj an (tu, T) K{t, T, uj) , K{t, T, lu) 
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sinh(^) 

where the two-point correlation function Gh{t,T) = ^^( Jh{t,x) J]^(0,0) )t can 
be computed on the lattice. Jh = <iij,x)VHq{T,x) is a local mesonic operator and 
r_H" = 7i,75 for vector {Vu) and pseudo-scalar {PS) channels, respectively. The 
temperature T is related to the Euclidean temporal extent aN^. by T = l/{aNr), where 
a is the lattice spacing. We measured charmonium correlation functions on very fine 
(a = O.Olfm) quenched lattices with a relatively large size of 128^ x 96, 128^ x 48, 
128^ X 32 and 128^ x 24 at 0.73 T^, 1.46 T^, 2.20 T^ and 2.93 T^, respectively. The lattice 
parameters and part of the results have been reported in Ref. [9]. 

2. Default model dependences of MEM results 

The Maximum Entropy Method (MEM) is used to extract spectral functions from 
correlators through Eq. ([1]) [10]. The input parameter of MEM is a default model 
(DM), which includes the knowledge about the spectral function, e.g. cr(a;) > 0. The 
output spectral functions from MEM are reliable only if they show a small dependence 
on the input DMs. Thus it is very important to study the default model depedence of 
the output spectral functions. 

From the left panel of Fig. [1] we see that the default model dependence of vector 
spectral functions at 0.73Tc is very small. The location of the ground state peak is very 
close to the physical J/ip mass. The width of the peak is too wide to be interpreted 
as the physical width of J/ip. We found no evidence for the existence of a transport 
peak at this temperature. In the right panel of Fig. [T]we show {G{tT) — GreciTT)) = 



^o((o)/to V|i @ 0.73 Tj 

DM1 
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Figure 1. Left: Default model dependences of vector spectral functions at 0.73 Tc- 
Right: differences of G{tT) from Grec{TT) as a function of tT. Grec{TT) is the 
reconstructed correlation function from the spectral function at 0.73Tc [9j. 
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J du)K{TT){a{uj,T)—a{uj,0.73Tc)) as a function of rT. The flatness of G—Grec at lAGT^ 
indicates some small changes in the bound states and the rising feature with increasing 
distances at higher temperatures indicates (o"(a;,T) — a{ijj,0.73Tc)) is negative in some 
low frequency region. To look into the detailed change of the spectral function from 
below to above Tc, it is crucial to investigate the spectral function itself. 

In Fig. [2] we show the default model dependence of vector spectral functions at 
T > Tc. We first vary the low frequency (transport) part of the spectral function in the 
default model, where plots Al) and A2) show the low frequency and high frequency parts 
of the spectral function, respectively. A small default model dependence is observed. We 
then vary the resonance part of the spectral function in the default model, where plots 
Bl) and B2) show the high frequency and low frequency parts of the spectral function, 
respectively. In particular, the first peak location in "DM 2" corresponds to the peak 
location of the spectral function at 0.73Tc. We observed that the ground state peak 
location in general shifts to higher energy region at lAQT^ and becomes fiat at higher 
temperatures. Transport peaks have relatively strong default model dependences as 
seen in plot B2). 

3. Conclusion 

We summarize current uncertainties of transport peaks in the left panel of Fig. H] and 
resulting heavy quark diffusion coefficients in the right panel. We found that 2ttTD 
at 1.46 Tc is close to unity and is slightly increasing with temperature. We also 
performed the default model dependence study in the PS channel and together with 
the uncertainties of vector spectral function are shown in Fig. |3l which suggests the 
dissociation of both J/tp and rjc at T ^ 1.46 Tc. 
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Figure 2. Default model dependences of vector spectral functions at T > Tc. The 
dotted and corresponding solid lines are for the different default models and output 
spectral functions, respectively. 
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Figure 3. Uncertainties of output spectral functions in PS (left) and Vu (right) 
channels at all available temperatures. The shaded areas are errors of output spectral 
functions from Jackknife and the solid lines inside the shaded areas are mean values 
of spectral functions. 
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Figure 4. Transport peaks at T > Tc and resulting diffusion coefhcients. 



